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Abstract 

At high energies, in particle-capture processes between ions and atoms, classical kinematic re- 
quirements show that generally double collision Thomas processes dominate. However, for certain 
mass-ratios these processes are kinematically forbidden. This paper explores the possibility of 
capture for such processes by triple or higher order collision processes. 
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I. INTRODUCTION 



It has long been known that at high energies, particle capture (also called particle or 



undamental three- 
J In 1927 L.H. 



mass transfer or exchange) is dominated by classical kinematics. As a 

body process it has been much studied in both experiment and theory 
n 

Thomas [5( did a theoretical analysis of the experiments of Rutherford, and especially G.H. 
Henderson js], in which energetic alpha particles captured an electron as they passed through 
several different media, emerging as He + ions. 

We denote the generic process 1 + (2, 3) — > (1, 2) + 3, and "high energy" means that the 
energy of the incident particle 1 is large compared with the binding energy of particle 2 to 
either the initial nucleus 3 or the final nucleus 1. Thus, without loss of significance we can 
consider the transferred particle 2 to be initially at rest with respect to the nucleus 3, and we 
define "capture" to mean that 2 emerges from the collision with zero velocity with respect 
to the scattered incident particle 1, i.e. they move off with identical speeds and directions. 
The nature of the interactions between the particles, generally Coulombic, is not important 
in our considerations, so we consider only hard contact collisions. 

Conservation of overall energy and momentum forbids the capture process to take place 
with only a single collision (except in the case where 1 and 3 have exactly the same mass 
and 1 collides head-on with 3). In the case Thomas considered, there are two collisions: In 
the first, particle 1 (the alpha particle) collides with the bound particle 2 (the electron). 
This collision brings the electron from rest to the speed needed for capture, but its direction 
is wrong. The electron then collides essentially elastically with the nucleus 3 to which it 
was bound, changing its direction to parallel that of the recoiling particle 1. The scattering 
angle of particle 1 is called the "Thomas angle" and is given by [1712/ mi) sin 60°. In the 
most studied case, protons on hydrogen, this is 0.472 mrad. We call this scenario "Thomas 
process A." 

Logically (but not necessarily physically) two other double collision processes are con- 
ceivable. In scenario "Thomas process B" the first collision between 1 and 2 brings particle 
2 to the necessary final speed. Particle 1 then collides with particle 3 and has its direction 
changed to enable capture of 2. In scenario "Thomas process C" the first collision is between 
particles 1 and 3. Recoiling particle 3 then collides with particle 2. These three possibilities 
are shown schematically in Fig. 1. Note that in all three cases, the kinematics are fixed 
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since the number of constraint equations (conservation of momentum and energy for each 
collision) equals the number of variables. However, depending upon the mass ratios, not 
all of these processes may be allowed. That is, conservation of energy and momentum may 
forbid one or more of these scenarios — or, in the cases studied below, all three may be 
forbidden. 

In 1987 the senior author of this paper devised a simple diagram to show how the available 
processes depend on the mass ratios |7| ■ This diagram, which has come to be known as the 
"Lieber diagram," is shown in Fig. 2. The symmetry of the diagram about the 45° line may 
be attributed to time-reversal symmetry. The curved boundaries are simple rectangular 
hyperbolas. Related diagrams have been given in 0, 

The quantum mechanical picture is more complicated [l-4|. Because of the uncertainty 
principle and the necessity of using wave packets, the momenta of the particles and their 
positions are spread out. At the high energies considered here, the Born approximation might 
be considered a good approximation. With Coulomb potentials it leads to a differential cross 
section with a large peak in the forward direction. The total cross section falls off as E~ 6 . 
However the second Born term dominates the first at sufficiently high energies because it 
falls off only as E~ 5 ^, as shown in the dissertation of Drisko [9]. This unique phenomenon 
occurs because of the Thomas double-collision process, which gives rise to a pole in the 
integrand of the second Born term when the propagator is on shell. The singularity gives 
rise to a second peak in the differential cross section at the Thomas angle which emerges 
from the background when the projectile energy is sufficiently large. The Thomas peak 
has been observed experimentally [lOj] in the pure three-body case of electron capture in 
proton-hydrogen scattering. An earlier atomic physics observation utilized capture of 
one electron from helium in proton helium scattering - the scenario of Rutherford, Henderson 
and Thomas - and so is somewhat less clear-cut because of the presence of two electrons. 

Even earlier there was an observation of the Thomas peak in a molecular collision, as 

n 

proposed by Bates [12J: a proton is captured from a methane molecule in the process p + 
CH4 — ¥ CH 3 + H2 '. The Thomas peak occurs at 46 degrees. 

The question to be explored in this paper is the behavior of collisions for which all three 
Thomas double-collision processes are kinematically forbidden. Is capture possible by more 
than two collisions? We answer this question in the affirmative. Atomic processes which 
are Thomas-forbidden are rare because of the particular masses of the electron and proton. 
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However, if we admit muons, the process: H+jj, — > e + (/ip), which describes the formation of 
muonic Hydrogen, can be shown to be Thomas-forbidden. Molecular processes are simpler 
to come by, e.g. Na + I 2 — > Nal + I. These and other processes are discussed below. 

Our study reveals interesting structure in the forbidden regions of the Lieber diagram. 
We have not determined the quantum mechanical behavior of the cross section for these 
processes, which would correspond to third and higher order terms in the Born series. In 
this situation there is no on-shell pole in the second Born term, but may be one in the third 
or higher term. It is not clear whether such a pole would lead to dominance of this term in 
the high energy cross section. 

In section II we briefly review the kinematics of the Thomas double collision processes. 
This will establish the notation and methods to be used for exploring the forbidden processes. 
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II. THOMAS DOUBLE-COLLISION KINEMATICS 



It is convenient to work with the ratios of the masses of particles 1 and 3 to the mass of 
particle 2: a = m\jmi and b = m^jm^- The mass ratios a and h which lead to scattering 
may be found simply by solving the equations of conservation of energy and momentum 
at each collision. It is also convenient to refer to the binary collisions using the number of 
the particle not involved in the collision, the "spectator particle." We will use parentheses 
to distinguish between particle numbers and collisions numbers. So, for example, (3) refers 
to a collision between particles 1 and 2. A capture process will be described by a string 
of numbers enclosed by parentheses corresponding to the collisions in the order they occur 
from left to right. Thus the three double collision processes we have considered, A, B, and 
C, can be referred to as (31), (32), and (21) respectively. A sketch of their trajectories is 
given in figure 1. 

We will now determine the mass ratios a and b for which these processes are allowed. 

First we look at process A, that is, (31). The initial velocities of particles 2 and 3 are 
zero and u will be the initial velocity of particle 1. After the first collision, (3), the velocities 
of 1 and 2 are u' and v. After the second collision, (1), the velocities of 2 and 3 are u' and 
w. The equations of conservation of momentum for the two collisions are: 

au = au' + v (1) 
v = u' + ow. (2) 

The equations of conservation of energy are: 

2 1 /2 , 1 2 /n\ 

-au = -au H — v (3) 
2 2 2 v ; 

\v 2 = \u' 2 + \bw\ (4) 

While solving these four equations is straightforward, we will use a technique which will 
make the calculations in Section III simpler. 

Without loss of generality we may consider the collisions to occur in a plane, the xy-plane. 
If we express the velocities of the three particles as complex numbers and write them as the 
components of a three dimensional, complex vector, as: (v\ x + ivi y , v 2x + iv%y, + 2t>3<J T 



then the velocities before and after a collision can be related by a 3 x 3 complex matrix 13 ]. 



If we denote the matrices associated with collisions (1), (2) and (3) by Si, S 2 , and 5*3, then 



FIG. 1: The three binary processes. (Bold line = particle 1, Thin line = particle 2, Dashed line = 
particle 3) 
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we can write: 



S^x) = S 1 (0) a ^ = 



1 + b 



1 + 6 
1 + bx b(l — x] 
1 — x b + x 



\ 



(5) 



s 2 ( x ) = s 2 (or^ = 



a + b 



a + bx 6(1 — x) 

a + b 
a(l — x) b + ax 



s 3 (x) = s 3 (or^ 



1 + a 



a + x 1 — x 
a(l — x) 1 + ax 
1 + a 



\ 



(6) 



(7) 



where a is the scattering angle of the collision and x = exp(ia). 

The matrix which transforms the initial velocities of the three particles to the final ve- 
locities is, for process A, Si(x 2 )S 3 (xi), where x\ = exp(iai) and x 2 = exp(ia 2 ) and a\ and 
a 2 are the scattering angles for the first and second collisions. The initial three velocities 
are: (u, 0, 0) T . By insisting that after the two collisions, 1 and 2 have the same velocity, we 
obtain an expression relating oi\ and a 2 : 

/A 



(l,-l,0)S 1 (x 2 )S 3 (x 1 ) 



u 





0. 



(8) 



Its solution is 



x 2 = 



[l + a + b)xi + ab 



(9) 

(ab-(l + a + b))/2ab. 



ab(l — xi) 

The condition that x 2 have unit absolute value implies that cos(a:i} 
We need cos(o;i) to be in the range [—1, 1] and this condition is satisfied when 3ab— a— b— 1 > 
0. Thus the allowed region for process A lies to the right of the rectangular hyperbola given by 
3ab — a — b — 1 = 0. The equation for the hyperbola can be written as (a — 1/3) (b — 1/3) = 4/9, 
so the asymptotes are at a = 1/3 and b = 1/3. 

For process B, or (32), the condition for capture is: 



{l,-l,0)S 2 (x 2 )S 3 (x 1 ) 



u 




\ / 



= 0. 



(10) 



FIG. 2: The Lieber diagram for double collision processes 




whose solution is: X2 = a(b — (1 + a + b)xi)/b(a + Xi). The absolute value of X2 is unity when 
cos(ai) = (a 2 + ab + a — b)/2ab. The cosine falls in the desired range when 3ab + a 2 + a — b > 
and b > a. Thus the allowed region for B lies between the line b = a and the curve 
3ab + a 2 + a — b = 0. The latter is again a rectangular hyperbola with an asymptote at 
a = 1/3 in the first quadrant. 

The process C, or (21), is similar to B. If B is viewed under time reversal and with 
particles 1 and 3 exchanged, it is equivalent to C. Thus the values of (a, b) for which C is 
allowed are simply those of B with the roles of a and b reversed. Thus C is allowed if (a, b) 
lies in the region between the the line b = a and the curve 3ab + b 2 + b — a = 0. The latter 
is a rectangular hyperbola with an asymptote at b = 1/3. 

The allowed regions for the three binary processes are shown in the "Lieber diagram," 
figure 2. 
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III. CAPTURE PROCESSES WITH THREE OR MORE COLLISIONS 



A. Preliminaries 

In order for a collision sequence to lead to capture, a few simple conditions must be met: 

1. Consecutive collisions must be distinct since two isolated particles cannot collide more 
than once. 

2. The first collision cannot be (1) because particle 1 initiates the capture process and the 
last collision must not be (3) because if 1 and 2 collide at the end of the process they cannot 
emerge with the same velocity. 

3. The times between collisions must be positive. 

Using these conditions, we prove the following useful lemma in Appendix A. 

Lemma 1 When three distinct collisions between three particles occur consecutively, a 
fourth collision must be of the same type as the second. Furthermore, sequences of the 
form (31323) in which the first, third and fifth collisions are of the same type, but the second 
and fourth are of a different type, cannot occur. 

Since capture begins with 2 and 3 having the same velocity and ends with 1 and 2 having 
the same velocity, the collision sequence can be considered to begin with the collision (1) and 
to end with the collision (3). If three distinct collisions were to occur consecutively between 
the virtual (1) and virtual (3), our lemma would be violated. Therefore a capture process 
must only consist of two types of collisions. There are only five such processes, namely: 
(32) n , (2)(32) n , (21) n , (21) n (2), and (31) n . This is fortuitous because the classical scattering 
matrices for such processes are computed simply by raising to a power the scattering matrix 
of the repeated pair. 

We will first consider (32) n and (2)(32) n , which have allowed regions which are within the 
forbidden region of the Lieber diagram. Because each collision in these processes involves 
particle 1, the whole problem can be done in the rest frame of 1. The scattering matrices are 
easily found for such a frame, but since particle 1 accelerates when it collides in an inertial 
frame and in the frame we are describing, it is at rest, this frame is noninertial. Thus energy 
and momentum will not be conserved in collisions (2) and (3) (the ones in which particle 
1 is involved). In this scheme, one simply applies the scattering matrices to the velocity 
vector and then subtracts from each of the three velocities, the velocity of particle 1. If we 
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have a matrix M and an initial velocity vector vl, the final velocity in the particle 1 frame, 
v} is given by: 



v f 



I - 



{A 

1 

W 
M 

1 

V) 



( 1 j Mv, 



(100) 



Mvi 



So to get the correct matrix one multiplies the matrix M by 



/ ^ 
-1 1 
-1 1 



V 



(11) 



(12) 



The benefit of working in the frame of particle 1 is that only two velocities need be dealt 
with and consequently the scattering matrices become two by two matrices, the first row and 
first column being irrelevant. The scattering matrices, 52 and S3, are in this representation 
given by: 



S2 



(13) 



1 b(y-l)/(a + b) \ S3= l x 

y ) ' \(x-l)/(a + l) 1 

One needs to know what the scattering angle is at each step to find the xs and ys. To 
that end, the following lemma, proved in Appendix B, is useful. 

Lemma 2 In a process of the form (2)(32) n or (32) n 7 all scattering angles, excluding the 
first, are equal. 
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B. Collisions of type (32) n 

The first step in analyzing (32) n or (2) (32)™ is to find the relation between the first and 
second scattering angles. Let the first and second scattering phasors be x = exp(ia) and 
y = exp(i(3). Suppose the particles 2 and 3 initially have unit velocity. Then following the 
first collision, the velocities of 2 and 3 are x and (a + x)/(a + 1). The second collision is 
between 3 and 1 and we should like to know its scattering angle, 0. If we exponentiate the 
equation derived above for scattering angle, we get: y = exp(i/3) = — exp(2iip) / exp(2i(p) . 
In the present case, exp(2i0) = [(a + x)/(a + l)]/[(a + x)/(a + 1)]* = x(a + x)/(ax + 1) and 
exp(iijj) = x, so that: 

y = - x (ax + l)/(a + x). (14) 

Now there is only one unknown variable in the problem, namely x. It is found by insisting 
that, in the particle-1 frame, 2 have zero velocity finally. To find 2's final velocity, we need 
to multiply the string of scattering matrices. Since only the first collision has a unique 
scattering angle, the latter 2n — 2 collision matrices may be expressed as a power to which 
the first two scattering matrices are appended: 

S = (S 2 S 3 y n - l) S { / /7r) si a/7r) = (S 2 S 3 ) n si a ^ y \ (15) 

To exponentiate S2S3, we need to diagonalize it. 

^ = .V + KV- 1)7(« + b){a+ 1) Ky - I)/ {a + b) 
y(y-l)/(a + l) y 



The trace and determinant of the matrix are given by: 

%-l) 2 



Trace(S 2 S 3 ) = 2y 



(17) 



2y(a + l)(a + b)_ 
det(S 2 S 3 ) = y 2 . (18) 

Since the determinant is the product of the eigenvalues, the eigenvalues may be expressed 
as ei = y/g and e 2 = y ■ g, where g is to be determined. Furthermore, since the trace is the 
sum of these, we have: 



2y l -{g + l/g) = 2y 



\ + Kv - if 



2y(a+l)(a + b) 



(19) 
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This equation suggests that we write g = e ld with 9 defined by: 

cos(0) = 1 + b ^~y 2 \ * (20) 
V ; 2y(a+l)(a + b) V ; 

It will now be shown that 9 is real and that the substitution is justified. From the identity 
(y - If/y = -4sin 2 (/3/2), we can write cos(#) = 1 - 26sin 2 (/3/2)/((a + l)(a + &)), so that: 



1 > 



(a + l)(a + b) > b 
6sin 2 (/3/2) 



> 



(a+l)(a + 6) ~ (a + l)(a + 6) 
X<1 26sin 2 (/3/2) ^ 



(a+l)(a + 6) 
-1 < cos(0) < 1. 

So 6* is real and the eigenvalues may be expressed as e± = y exp(±i9) . S may be decomposed 
into idempotent matrices, P+ and P_ as S = e + P + + e_P_, if P + = (S — e_J)/ (e+ — e_) and 
P_ = (5-e + /)/(c_-e+). 

We are now in a position to state the condition for capture. The final velocity of particle 
2 is 0, so the condition is: 



1 J (S 2 S 3 ) n S : 
Inserting the eigenvalue expansion, we have: 



(a-/3)/7r 
3 



0. 



(21) 



0=10 



1 



S 6—1 S — tjl-I 

e™ h e n — 



e + -e_ 



e_ — e_ 



' y " _1 ? [e me (5 - ye- 19 /) - e- mB {S - ye ie I)] S^ )/7r ' ' 



1 



e ie - e 



,n—l r 



1 



sin(0) 



sin(n0)sf /7r) sf /,r) -|/sin((n - 1)0)/ 



5. 



(a-/3)/vr 



When the appropriate matrix elements of s!f ^'S^ and S 3 
Eq. ([21]), we obtain: 

sin((n — 1)9) ^ 26(acoso; + l) 



(a-/3)/7r 



(22) 
(23) 
(24) 



are inserted into 



sunt 



(a + 6)(a+l) 



(25) 
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It will be useful to eliminate the alpha dependence from this equation. First we invert 
Eq. (TH|) . As it is quadratic in x, there are two solutions: 



x 



-(1 + y)± v / (l + y) 2 -4a 2 ? /l /(2a) 



cos a 



COS' 



! (/3/2) ± sin(/3/2) Va 2 - cos 2 (/3/2) /a. 



(26) 
(27) 



If either solution is inserted into Eq. (I2~5"j) . and Eq. ( 12U|) is applied, we obtain the following 
relation: 

cos 2 (fl/2) q(l + 6) 

sin 2 (^) (a + 6)' 1 J 

Now for given a and 6, there may be several solutions to Eq. fl28l) . We can determine their 
legitimacy by looking at the signs of the time intervals between collisions. We construct a 
formula for the ratios of successive times. Consider the initial positions and velocities of 
the particles. If the process occurs in the x-y plane, with 2 at the origin and l's velocity in 
the positive x direction, and if the polar angle of 3's position vector is <fi and the distance 
between 2 and 3 is unity, then the initial position vector is ^ j ■ The position vector 
at the moment of the (2j) th collision is: 



r 2j 



1 x + t 2 s^st h [ 1 I + 



(29) 



■t 2 j-\ o 3 b 2 



At the moment of the (2j + 1) it is: 




r 2 i+i 




, Q /3/n Q a/ir 




(30) 



l 2 j I o 2 D 3 



The position vectors of relation (|29|) must have their second component zero because 
particle 3 is at the origin at the moment of an even numbered collision. Similarly, the 
position vectors of relation (|30|) must have their first component zero because particle two 
is at the origin at the moment of an odd numbered collision. From these requirements, we 
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can write equations from which to determine the times between collisions. Their basic form 
is: 

( 1 ) r 2j = and ( 1 ) r 2j+1 = 0. (31) 

First we will find the ratio between an even numbered time and the preceding odd numbered 
time. If we subtract the equation ^ 1 j r 2j-i = from the equation ^ 1 j r 2j+i = 0> 
applying Eq. f l30|) . we find that: 



o = V-i (io) (s^s^y'st^ 1 ' 



(32) 



+t 2j [io)[s^s^Yst" )/n [ I 



If we apply the identity ( 1 j S^ w v — y ( 1 j v, which holds for any velocity vector v, 
and carry out the exponentiations, we find: 



1 



sin((j - l)6)S^ M St M - y sin(C7 - 2)0)S< 



l2j 



1 ) sin(ms¥ M St M -ysin((j -l)6)St" )/n 

1 



(33) 



The ratio of the two matrix elements of the numerator and denominator can be extracted 
from Eq. (122"]) : 



( i o ) sf^st M I j 



(i o)s : 



3 



sin((n - 1)0) 
sin(n0) 



If this ratio is applied to Eq. (I3"3"j) . it can be shown that: 
tjy sin((n- j + l)0) 



<<2j-l 



sin((n — j)9) 



j = l,2,...(n-l) 



(34) 



(35) 



It remains to find the ratio of an odd numbered time to the preceding even numbered 
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time. We subtract ^ 1 J *2j = from ^ 1 J r 2j+2 = 0, to find: 



o = t 2j [ o i ) (s^s^Y sf* ' 1 



(36) 



+t 2j+1 ( 1 ) (s^s, 



i+i / / f 



Following the same procedure as above, we use the fact that ^ 1 j is a left eigenvector of 
and that S3S2 and S2S3 have the same eigenvalues to derive: 



1 



sin((j - 1)9)S^S^S^ -ysm((j - 2)9)S^ 



l 2j+l 
t 2 j 



(37) 



1 



sm(je)S^StS^ -yma(U - 1)9)S^' 



As before, we first evaluate the ratio of the matrix elements, p. If we perform the matrix 
multiplications and apply Eqs. (|14|) and (|20p . we find that: 



P 



(0 1 ) sf^s^s*? 1 ^ r 

(oi)^^j 



2cos# - 1 



2 (a cosa + 1) 
a 2 + 2a cos a + 1 



(3? 



Next we solve Eq. ( 125]) for 2(a cos a + 1) in terms of 9, a, and 6, and use the result to simplify 
P ■ 



p = y 



2cos# - 1 



sin(n#) — sin((n — 1)9) 



(39) 



(40) 



[a(b + I) /{a + b)} sm(n6) - sin((n - 1)9) 
Finally, if we eliminate a(b + l)/(a + 6) using Eq. (1281) . we obtain: 

p = ycos((n - 3/2)9)/ cos((n - 1/2)0). 

And when the expression for p is used in Eq. (137)) . it can be shown that: 

t 2j+1 _ COs((n-j + 1/2)0) 

t 2 , cos((n-j- 1/2)0)' 3 ^ 

As was mentioned earlier, for a given a and 6, there may be several values of 9 which 
satisfy Eq. (|28|) but at most one corresponds to a process with positive times between 



(41) 
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collisions. The time between the first two collisions, t\ can be guaranteed to be positive 
because after 1 collides with 2, a collision with 3 is possible if 3 has the right orientation. 
But from relations ( 13 5 p and (141 p . we see that for all the times to be positive, we need 
sin((n- j + l)0)/sin((n- j)0) < and cos((n- j + 1/2)0)/ cos((n -j - 1/2)0) < for all 
1 < j < n - 1. 

The solution set for sin(/c0)/ sin ((A; — 1)0) < is 

fc-i 

\J(ln/k,ln/(k-l)) (42) 
i=i 

and the solution set for cos((A; + 1/2)0)/ cos((A; - 1/2)0) < is 

k 

\J((2l - l)n/{2k + 1), {21 - l)7r/(2fc - 1)). (43) 
i=i 

In (|42p and fj43|) . k = n — j and j varies from 1 to n — 1 so fc can take any value between 
1 and n — 1. It is not hard to show that the intersection of the sets (142]) having A; between 
1 and rz — 1 is ((n — 2)ir/(n — 1), 7r). The intersection of this interval with all the sets (I43p 
having A; between 1 and n — 1 is ((2n — (2n — 1), 7r). Thus a necessary condition for all 
positive time intervals is that be in the range ((2n — 3)7r/(2n — 1), 7r). 

We have seen that is real if j3 is. From the form of Eq. ( I14p . the reality of a guarantees 
the reality of /3. 

All the dynamical variables of the problem are found by multiplying the original velocity 
vector by the 2x2 matrices S2 and S3. So a process is classically allowed if a G 5i and 
{2n — 3)7r/(2n — 1) < < 7T. Eq. ( 126|) may be rewritten as x = A /^exp(i0), if cos0 = 
— (1 + y) / (2a,y/y) = — cos(/3/2)/a. In this form, it is evident that if (3 e 3?, then a G 3?, iff 
cos 2 (/3/2) < a 2 . This inequality may be recast in terms of sines as: 

1 - a 2 < sin 2 (/3/2) < 1 ^ ~ " } < sin 2 (0/2) < 7 -. (44) 

~ v / 1 (a + b) ~ w 7_ (a + 6)(a+l) v ' 

The lower bound on sin 2 (0/2) poses no constraint, it being guaranteed by Eq. (T28l) . so we 
must solve that equation in the domain (2n — 3)ir/ (2n — 1) < < 7r, subject to the constraint 

sin 2 (0/2) < 6/((a + 6)(a+ 1)). (45) 



Let Q(0) be defined by Q{9) = cos 2 (0/2)/ sin 2 (n0). Then Eq. can be written as 

jl + j 
a + b 
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» = W^'l. («) 



This equation can be used to assign a value to at every point in the ab plane at which 
Q^ 1 exists. On the interval ((2n — 3)7r/(2n — l),7r), Q{9) decreases monotonically from oo 
to 1 /An 2 . Thus, we may only choose a and b for which 

> 1 /An 2 . (47) 



a (! + b ) , ■■ / ■ . 2 



a + 6 

The points (a, 6) which satisfy the inequality a(l + b)/(a + b) > Qq lie between the branches 
of a hyperbola which passes through the origin and has asymptotes at a = 1/Qq and b = 
l/Qo — 1. Because Q decreases monotonically on the interval of interest, if 6i, 02 and 83 
satisfy (2n - 3)n/(2n - 1) < 8 X < 2 < 8 3 < tt, it follows that if 8 X < 8{a,b 1 ) < 8 2 and 
82 < 8(a, 62) < #3, then 6 2 > In other words, the region in which 8\ < 8 < 8 2 lies under 
the region in which 82 < 8 < 83. 

We can find the allowed region for this process by looking for the curves along which the 
left and right hand sides of the inequality f|4"5]) are equal. Then we can test to see on which 
side of the curves the inequality is met. Finally we must check the inequality (|47|) to make 
sure that 8 actually exists in the proposed region. 

We write H5] as an equation: 

sin 2 (fl/2) = 6/((o + 6)(a + l)). (48) 
and combine it with Eq. (I2"8"|) to obtain: 



> 1 + a + b 



tan 2 (n#) = (49) 
ab 



Also, Eq. fl48|) can be modified to yield: 



tan 2 (0/2) = -j— -. (50) 

w ' a(l + a + b) v ' 



If we multiply Eqs. fl49|) and fl50|) together and solve for a, we get two solutions. The two 
solutions for a along with the corresponding expressions for b are as follows: 

a = cot(0/2)cot(n0) (51) 
b = cot(n0) cot((n - 1/2)0), if (52) 
(2n - 2)7r/(2n - 1) < < (2n - l)vr/(2n) 



and 



a = -cat(0/2)cat(n0) (53) 
b = cot(n0) cot((n + 1/2)0), if (54) 
(2n - l)vr/(2n) < < (2n)7r/(2n + 1). 
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The limits on 9 are determined by finding the values of 9 which make a and b both positive. 
For convenience, let the former be denoted by C 2n -\ and the latter by C 2n . If we let 
9 X = (2n - 2)vr/(2n - 1), 9 2 = (2n - l)?r/2n and 9 3 = 2nir/(2n + 1), then it follows from 
the above discussion that C 2n _\ lies beneath C 2n . By plugging 9 — 9\ into f[5Tj) . we find 
that C 2n -i has an asymptote at a = tan 2 (7r/(4n — 2)) and by plugging 9 = 9 3 into (IB^|) . 
we find that C*2 n has an asymptote at a = tan 2 (7r/(4n + 2)). To test whether the process 
is allowed between these two or outside them, we study a test point. If b = and a > 0, 
then 9 exists and is given by Q~ l (l). Since (2n — 3)7r/(2n — 1) < 9 < tt, it follows that 
sin 2 (6*/2) > cos 2 (7r/2n) but from inequality fH5"|) we need sin 2 (6'/2) < at such a point. 
Thus the allowed region for the process is between the curves C 2n -i and C 2n - 
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C. Collisions of type (2)(32) 



The next process to consider is (2) (32)". The relation between the first and following 
scattering angles is different for that problem than for (32) n . In the above derivation of the 
relation, only the effects of the first two collisions are considered so the equation for the 
(2)(32) n problem is obtained by exchanging the masses of 2 and 3 in Eq. (|T4|) : 

y = —x(ax + b)/(a + bx) . (55) 

The condition for capture is: 



i o ) (s 2 s 3 ) n s^ ^ I j 



1 



wn(n8)S% /ir Sl /v S% /ir -ysm{{n - l)6)S, 



a/ir 



l 



(56) 



(57) 



If we evaluate the matrix elements with the aid of Eq.( 1201) . we find that: 

sin(n60 a 2 + lab cos a + b 2 



(58) 



sin((n + 1)0) a 2 - b 2 

A formula is needed for a given fl. This is obtained by exchanging the masses of 2 and 3 in 
Eq. ([26]): 



cos a 



-6cos 2 (/3/2) ± Va 2 -b 2 cos 2 OS/2) sin(/3/2) /a. 
This formula, along with Eq. fTSOj) allows us to rewrite Eq. ( 158]) as: 

cos 2 ((n + 1/2)6) 6(1 + a) 



(59) 



(60) 



cos 2 (fl/2) a(l + b) 

As with (32) n , we next determine which values of 9 can be ruled out because of negative 
time intervals between collisions. Since 3 begins at the origin, will now denote the initial 
polar angle of 2. The position vector at the moment of the (2j) th collision is: 



J4> 



r 2j 



+ hs: 



(61) 



•hj-x [s^S- 
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At the moment of the (2j + 1) it is: 



r 2j+i 



+ hS. 



+ t 2 Sg /v S% /ir 



+ 



(62) 



■ ■ t 2j [ S 3 S 2 



As before, we derive the time ratios from the constraint that the two particles involved in a 
collision have the same position. The equations expressing that constraint are: 



( 1 ) r 2j = and ( 1 ) r 2j+1 = 0. 



(63) 



Now we can use Eq. (1621 to find the ratio of an even numbered time to the preceding odd 
numbered one: 



*2j 



1 



sin((j - l)9)StsT -ysin((j - 2)9)S { 2 



^2j-l 



1 



sm(je)S^ST-ysm((j-l)e)S ( 2 



The relevant matrix element ratio is: 



(64) 



1 ) S§ /7r S^ 



P 



1 s. 



-,(a-/3)/vr 



2(6 + a cos a) 
(a + l)(a + 6) 



(65) 



With the help of Eqs. ( 158|) and ( |60i) . we may eliminate a, b, and a from this expression: 
p = ycos((n — 1/2)9)/ cos((n + 1/2)9). When this is used in Eq. fl64"]) . we obtain: 



cos((n- j + 3/2)0) 
t 2 j-i ~cos((n- j + 1/2)0)' 



'2; 



j = 1,2, • • -n. 



(66) 



From Eq. (IBTj) . the ratio of an odd numbered time to the preceding even one is given by: 



y(lo)(s^S^y 1 S : 



l 2j+l 
t 2 j 



i o ] [S^s^Y 1 1 



(67) 
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Invoking Eq. (|56|) . we can remove the a, 6, and a dependence: 

t _ s ia{ (n -j + m . = 12 ..(„_!,. (6g) 

The solution set of cos((A; + 3/2)0)/ cos( (Jfe + 1/2)0) < -1 is 

fc+i 

lj((2Z - 1)tt/(2A; + 3), (2/ - 1)tt/(2A; + 1)). (69) 
and the solution set of sin((& + 1)9)/ sm(k9) < —1 is 

k 

\J(h/(k + l),h/k). (70) 
i=i 

In (|70p . k = n — j and j varies from 1 to n — 1 so k can take any value between 1 and 
n — 1. It is not hard to show that the intersection of the sets (I70p having k between 1 and 
n — 1 is ((n — 1)71/11, tt). In ( )69|) . j varies from 1 to n so k takes values between and n — 1. 
The intersection of ((n — l)7r/n,7r) with all the sets (169]) having /c between and rz — 1 is 
((2n — l)ix/{2n + 1), it). Thus a necessary condition for all positive time intervals is that 9 
be in the range ((2n — l)7r/(2n + 1), tt). 

Eq. ( 1551) may be inverted as x = ^expicj), if cos((/>) = —6(1 + y)/{2a^/y) = 
-6cos(/3/2)/o. We have G 3? if < cos 2 (/3/2) < a 2 /b 2 . Equivalently, we may write 
the condition for capture as: 

1 - a 2 /b 2 < sin 2 (/3/2) < 1 & §- - a \ < sin 2 (0/2) < -. -. (71) 

' v 6(o + l) ' ~ (a + b)(a+l) v ' 

The left inequality follows from Eq. ( 160]) so any solution to Eq. ( 160]) in the interval (2n — 
l)n/(2n + 1) < 9 < tc satisfying 

sin 2 (0/2) < 6/[(o + 6)(o+l)] (72) 

corresponds to allowed capture. Let Q be defined by Q{9) = cos 2 ((ri + 1/2)9)/ cos 2 (6 l /2). On 
the interval under consideration, Q increases monotonically in 9, from at the left endpoint 
to (2n + l) 2 at the right endpoint. From Eq. (1501) . we can show that 6 = aQ/(l + a — aQ). 
When this is used in the inequality we have: 

Q 



(l + Q) + o(l-Q) 



> (o + l)sin 2 (0/2). (73) 
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Taken as an equation, the two solutions of (173]) for a are — cot (9/2) cot (n0) and 
cot(6 l /2) cot((n + 1)0). The corresponding expressions for h and the values of which make 
a and b positive are as follows: 



a = - cot (0/2) cot (n0) (74) 
b = cot((n + 1/2)0) cot(n0), if (75) 
(2n - l)7r/(2n) < < 2n7r/(2n + 1) 



and 



a = cot (0/2) cot ((n + 1)0) (76) 
6 = co t((n + 1/2)0) cot((n + 1)0), if (77) 
2mr/(2n + 1) < < (2n + l)n/(2n + 2). 

The first curve is the same as C*2n and the second is the same as Cm+\- Thus, the boundary 
curves for 2(32) n are the same as those for (32) n . The functions a\ and <22 are found by 
combining Eq. f )60|) with the Eq. ( ITT]) . This has been done by eliminating b from Eq. (1711) . 

But a more useful form for the curves is obtained by eliminating 0. If Eq. (ITT]) is taken 
as an equality, then exp(i0/2) may be written as 

exp(z0/2) = -. (78) 

A/(a + 0)(1 + a) 

From Eq. ([60]), cos 2 ((n + 1/2)0) = b(a + 1 + fc)/((a + 6)(1 + b)) and 

exp(*(n + 1/2)0) = -. (79) 

y/(a + b)(l + b) 

The sign of sin((n + 1/2)0) depends on the parity of n. For even n, sin((n+ 1/2)0) is negative 
if (2n-l)7r/(2n) < < 2wr/(2n + l) and positive if 2mr/(2n + 1) < < (2n + l)vr/(2n + 2), 
while for odd n, the opposite is true. Combining Eqs. ( 178]) and ( 179]) gives for Ci n : 



2n+l 



and for C< 



2n+l- 



v/a(a + 1 + 6) +i v / &\ ^&(a + 1 + 6) - (-l) n iVa 

A/a(a + 1 + 6) - iVfe y ^(a + l + b) + (-l) n i^/a 

, ^\ 2n+l , 

v/a(a + 1 + b) +iVb \ ^Jb(a + 1 + b) + (-l) n i^/a 

y/a(a + 1 + 6) - iVb J ^b(a + 1 + b) - (-l) n i^/a 
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(80) 



(81) 



The first four boundary curves are given by: 



= b -a, (82) 

= 3ab-b + a 2 + a, (83) 

= 3ab 2 - b 2 + 2ba 2 + 6ab - a 3 - a 2 , (84) 

= 5a V - 10ab 2 + b 2 + 66a 3 - 6ab + a 4 + 2a 3 + a 2 . (85) 

These curves are, respectively, linear, quadratic, cubic, and quartic. Because the (32)™ 
and (2) (32)™ processes share boundary curves, the allowed regions are contiguous and non- 
overlapping. Furthermore, the asymptotes, having the basic form a = tan 2 (7r/(2(2m + 1))), 
approach a = in the limit as n approaches oo, and so the entire a — b plane above the 
line b = a is covered by allowed regions for processes of these types. The allowed regions for 
(32), (232), and (3232) are given in figure 3. 

FIG. 3: The pairs of curves which, going from right to left, give the allowed regions for (32), (232), 
(3232) and higher order processes. 




a 
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Starting at the a-axis and going counterclockwise, the figure shows curves C\, C2, C3 and 
C4. In the figure we can see that the pair C\ and C2 as well as the pair C3 and C4 have 
the same slope at the origin. Also C2 and C3 have the same asymptote. It is not difficult 
to show from Eqs. [SO] and EH that, i n general, C<m+\ and C*2n+2 have the same slope at the 
origin and C 2n and C 2n +i share an asymptote. 

A couple of examples of (2)(32) n type processes were mentioned in the introduction. 
These and one other example are given in more detail now. The capture of a proton from 
protonium by a muon is a (232) process: 

/I + (pp) -> (pp) + p. (86) 

The muon is particle 1, the proton is particle 2 and the antiproton is particle 3. The 
scattering angle can be shown to be 2.72 rad. The process Na + 12—)- Nal + I is a molecular 
example of (232) which has scattering angle 1.09 rad. The capture of a proton from muonic 
hydrogen by an electron: e+ (fip) — » H+ \i is an example of a higher order process, (2) (32) 10 , 
with scattering angle 1.09 rad (the same as the previous angle only to three significant 
figures). The time reversed process, H + //—)■ e + (up) is more interesting because it is a 
process in which muonic hydrogen is created. It is a (21) 10 (2) process with scattering angle 
59 mrad. 
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D. Collisions of type (21) n or (21) n (2) 

We know that at the beginning of a capture process, 2 and 3 have zero relative velocity 
and at the end, 1 and 2 have zero relative velocity. If we examine a given process under 
time reversal, we see that, initially, 1 and 2 have zero relative velocity and, finally, 2 and 3 
have zero velocity. If the particle labels "1" and "3" are exchanged, one has a new capture 
process. If this time reversal and label exchange is done to the process (32) n , what results is 
(21) n and when it is done to (2)(32) n , what results is (21) n (2). This means that if (32)" or 
(2)(32) ra is allowed for a point (a, b) = (a , b ), then (21)™ or (21) n (2), resp., is allowed for the 
point (a, b) = (b ,a ). Therefore the allowed regions for these processes may be generated 
by reflecting the previously determined regions about the line b = a. 



25 



E. Collisions of type (31) n 

The (31) n problem is much simpler than the processes previously considered because, as 
will be shown, each collision has the same scattering angle. Because of this there will be 
only two variables to solve for as opposed to the three in the above problem. Since 2 is 
involved in each collision, we work in its frame. The scattering matrices are: 

a=| lt(l - 1)/(1+6) V*=f v °l m 

x J \a{y- l)/(l + a) 1 

The condition which determines the relation between the first two scattering angles is that 
the second collision (between 3 and 2) reverse particle l's velocity so that it can collide with 2. 
Or in the case of (31), l's velocity is not reversed but made zero. In the frame of 2, the initial 
velocity is ^ 1 j , the velocity following the first collision is: ^ x a(x — 1)/(1 + a) ^ , 
and the final velocity is: 



(88) 



x + ab{x - l)(y - 1)/(1 + a) /{I + b) 
ay(x — 1)/(1 + a) 

Before the second collision, 1 has velocity x and following that collision, it has velocity 
x + ab(x — l)(y — l)/(l + a)/ (1 + b) so if l's velocity is to be reversed or brought to zero, x and 
(x — l){y — 1) must have opposite phase. The phase of the former is a and that of the latter 
is 0/2 + a/2 + it, so the first and second scattering angles are equal. This is demonstrated 
for the case n = 1 by Bittensky, through a geometrical argument jsj. We can therefore find 
the total scattering matrix simply by exponentiating M = S^ n . M, given by: 

M= I x + ~ l)7(a + !)/(& + 1) K* ~ 1)/(1 + b) . (gg) 
ax(x — 1)/(1 + a) x 

has determinant x 2 and trace 2x + ab(x — l) 2 /(a + + 1). As before, we express the 
eigenvalues as e + = xg and e_ = x/g, with g = exp(i9) so that 

06 f (x — 1) 2 \ . . 

cos (6) = 1 + — r ^ . 90 

V ; (o + l)(6 + l) \ 2x ) y ! 

The projection operators are: 



26 



Using these expressions and the above for e + and e_, we find that 



M n = e + P+ + e_P_ (93) 

™Tl— 1 

= — — (sin(n0)M-xsin((n- 1)6)1) (94) 

S11T ( u J 

The process is completed when 1 has zero velocity so the condition for capture is: 

( 1 ) M n ( * \ = =4* (95) 

xsm{(n-l)6)/sm(n6) = (lO^M^M (96) 

= x + ab(x-l) 2 /(a + l)/(b+l). (97) 

When this equation is combined with the definition of cos(6 l ), the simple result: 
sin((n + 1)9) = sm(n9), emerges. Aside from the trivial solution 9 = 0, which corresponds 



to no collisions at all, there are n solutions: 



9 = (2m + l)n/(2n + 1), m = 0, 1, 2, . . . , n (98) 

The allowed region is now readily found by insisting that a be real, which is to say, 
— 1 < cos(a) < 1. If we solve Eq. fl9U|) for cos(a) and insist that it be real and between -1 
and 1, we are left with: 



1 + a + b < abcot 2 (9/2). (99) 

From this equation, it follows that 9 cannot be it for that would yield the condition: 
1 + a + b < 0. However, the remaining n possible choices for 9 yield allowed regions which 
are bounded by the hyperbola: 1 + a + b = cot 2 (9/2). The hyperbola is symmetric about 
a = b and has asymptotes at a = tan 2 (6*/2) and b = tan 2 (6*/2). For example, if n=l, we 
are dealing with (32) which is the symmetric curve of the Lieber diagram. In that case, 
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6 = 7r/3 and our analysis gives the allowed region described by: 3ab — a — b — 1 > 0. To 
ascertain which 9 is the physical one, we insist that the time duration between collisions are 
all positive. The formulas for the time ratios here are: 

W*2i-i = " cos((j - 1/2)6)/ cos(0' + 1/2)0) and (100) 
ty+i/ty = -sm(j0)/sin((j + l)0). (101) 

These must be nonnegative for j ; < n — 1, which requirement holds only for (n — l)7r/n < 
6 < 7i. The only suitable from Eq. fl98|) is (2n — l)7r/ (2n + 1) and so the allowed region for 
(31) n is given by afrcot 2 ((2n - l)7r/(2(2n + 1))) - o- 6- l>0. That is, it lies to the right 
of the right branch of the hyperbola given by: 

a&cot 2 ( (2n-lM _ _ 6 _ = ( 2) 

V2(2n + l)y V ; 
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IV. CONCLUSIONS 



We have shown that there exist regions in the Lieber diagram for which the processes 
(2) (32)", (32)™, (21)™, (21)™(2) and (31)™ can occur classically and that, except for the last 
type, these regions are non-overlapping. If mass ratios are chosen in figure 3 which lie in the 
forbidden region, there is a unique process whereby particle 1 can capture particle 2. (The 
forbidden region excludes processes of the fifth type.) 



Appendix A: Proof of Lemma 1 



We will first demonstrate the impossibility of the sequence (2312). In what follows, the 
respective velocities of 1, 2 and 3 will be given in column vector form with v representing 
the three velocities and f representing the three positions. Let the velocities and positions 
for 1, 2 and 3 immediately after the first collision be 



I o N 

v 

y we ia j 



and 



l o n 



V / 



(Al) 



The quantities v, w, and r are positive reals and a is a real angle. We are working in the 
reference frame of particle 1. 

In sequences of the form (232), the scattering angle of the second collision is given by 
9 = 7i + 2a and in sequences of the form (231), the scattering angle is given by 9 = 
a + sin -1 ((sin a) /a). 

The time until the next collision is given by t± — r/v. At the end of this interval, the 
positions are: 



r = 



( o n 





(A2) 



rwe %a j v 

Let the scattering angle of the second collision be 9. Then after that collision the velocities 
are given by 

( v(l-e w )/(l + a) N 
v(l + ae ie )/(l + a) ■ (A3) 



v = 



\ 



we 



I 
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or in the reference frame of 2: 



-ve 




\ 



(A4) 



\ywe ia -v(l + ae id )/(l + a) J 

The time until the third collision is given by 

rwe ta /v 



to = 



(AS) 



[1 + ae ie )v/(l + a) - toe"' 

In order for t 2 to be real and positive, (1 + ae l6 )v/ (1 + a) must have the same phase as and 
be greater in magnitude than we ia . Assuming this to be the case, t 2 can be expressed as: 



rw/v 



(y/l + 2acos{6) + a 2 )v/{l + a) - w 
At the end of this interval the positions are 



t -ve i9 t 2 ^ 



(A6) 



r = 



V 



o 





(A7) 



/ 



Let the third collision have scattering angle <fi. Then after the third collision the velocities 
are given by 

( -ve ie ^ 

v= b(l-e^)(we ia (l + a)-v(l + ae ie ))/(l + a)(l + b) ■ (A8) 
v (b + e i(t> )(we ia (l + a) -v(l + ae ie ))/(l + a)(l + b) 

The time until the fourth collision is given by 

ve i9 t 2 



U = 



(A9) 



(6 + e^)(u»e ia (l + a) - v(l + ae w ))/(l + a)(l + b) +ve id ' 

In order for t 3 to be positive, (b + e i(f) )(we ia (l + a) - v(l + ae ie ))/{l + a)(l + b) must 
have the same phase as and be greater in magnitude than —ve ld . The magnitude of 
the former is \Jb 2 + 2b cos + l(vy/ a 2 + 2a cos + 1/(1 + a) — w)/(l + b). The factor 
vy/ a 2 + 2a cos 9 + 1/ (1 + a) — w has a value less than v and the factor \/b 2 + 26 cos + 1/ (1 + 
6) is less than 1 so the condition for the third time interval to be positive cannot be met. 
Therefore the proposed collision sequence is impossible. 
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We will next demonstrate the impossibility of processes of the form (32313). 

Immediately after the first collision, in the frame of 1, 2 is moving away from 1 and 3 is 
moving toward 1. Let the velocity of 3 be w and the velocity of 2 be vexp(ia). We view 
the next collision, (2), in the frame of particle 1. After the collision (2), 2 needs to have its 
velocity shifted by a phase of 7r so that it changes from an outgoing particle to an ingoing 
particle. The collision matrix for collision (2) is: 



(AlO) 



1 ^sin(0/2)exp[»(0 + 7r)/2] 
exp(i#) 

Before the second collision, the velocities of 2 and 3 are (v exp(ia), w) T . If we apply the 
above matrix to this vector, we find the velocities after the second collision: 

(vexp(ia) + 2bwsm(6/2) exp[i(9 + vr)/2]/(a + b), wexp(i9)) T . (All) 

As mentioned, we need to have the velocity of 2 shifted by a phase of tt. This requires 
that the term 26-u; sin(#/2) exp[i(9 + 7r)/2]/(a + b) have opposite phase from v exp(ia) and 
be larger in magnitude than v. For the phase of the correction term to be right, we need 
9 = 2a + 7T. Plugging this into Eq. IA11I gives the velocities as: 

((t> — 26u;cosa;/(a + b)) exp(za), — w exp(2ia)) T (A12) 

and for the velocity of particle 2 to be turned around, we need: 

2bw cos a /a \ 

v < — . A13 

a + b v ' 

We will analyze the third collision by working in the frame of 1 before the collision and 
by working in the frame of 2 after the collision. That way, before the collision, 3 is moving 
away from 1 and 2 is moving toward 1 and after the collision, 1 is moving away from 2 and 
3 is moving toward 2. At the moment of the third collision, 2 replaces 1 at the origin. Thus 
before the collision, 3 is moving away from the origin and after the collision it is moving 
toward the origin. So we need the collision to change the phase of particle 3 by n. The 
matrix for collision (3) in which we work in the frame of 1 before the collision and the frame 
of 2 after the collision is: 



— exp(i(p) 
-(1 +aexp(i0))/(l + a) 1 
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(A14) 



Before the the third collision, the velocities of 2 and 3 are : 

exp(ia)((v — 2bw cos a/ (a + &)), — w exp(ia)) T (A15) 
which may also be written as: 

— exp(ia)(V, wexp(ia)) T (A16) 

where v' is a positive number. From the matrix IA14| the velocity of 3 after the collision is 
— exp(ia)(w exp(ia) — (1 + aexp(i0))i// (1 + a)) the phase of which must differ from that of 
wexp(ia) by ir. This requires that (1 + aexp(i0))t>'/(l + a) must have the same phase as 
and be greater in magnitude than wexp(ia). For the latter to have the same phase as the 
former, we need 

(j) = a + sin 1 ( — — J . (A17) 

If the above value for is used, the magnitude of (1 + a exp(i(p))v' / (1 + a) can be shown to 
be (cos a + \J a 2 — sin 2 a)v '/ (1 + a). Since this must be greater than w, the condition for 
the fourth collision is: 



2bw cos a 

; v 

a + o 



(cos a + v a 2 — sin 2 a) . , , „. 

> w (A18) 

1 + a 



After the third collision, in the reference frame of 2, the velocities of 1 and 3 are: 

— exp(2ia)(— v'{y/ a 2 — sin 2 a + i sina)/a, w — (cos a + V a 2 — sin 2 a)v ' / (1 + a)) T . (A19) 

or 

- exp{2ia)(-v' exp{i(<j) - a)), -w') T (A20) 

where w' is a positive number. We can analyze the fourth collision in the reference frame 
of particle 2. Before the fourth collision, 1 is moving away from 2 and after, 1 is moving 
toward 2. Thus the fourth collision must change the phase of the velocity of 1 by 7r. The 
matrix which mediates this collision is: 



1 26exp[i(V> + 7r)/2]sin(^/2)/(l + 6) 
exp(i^) 



(A21) 



Before the fourth collision, the phase of the velocity of 1 minus the phase of the velocity 
of 3 is — a. Thus the scattering angle of the fourth collision must be 2(0 — a) + 7r and the 
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matrix may be rewritten as: 



1 — 2bexp[i((f) — a)] cos(0 — a)/(l + b) 
-exp(2i(0-a)) 

The velocity of 1 after the fourth collision may be written as: 

— exp(2ia)(— v' exp(i((j) — a)) + 2bw' exp[i(0 — a)] cos(0 — «)/(! + 6)) 



(A22) 



(A23) 



or exp(i(0 + — 2 frit/ cos(0 — a)/(l + 6)) whereas before it was exp(i(0 + a))i/. Thus 

the condition for the fifth collision is: 



2Wcos(0-a)/(l + 6) > u' 



(A24) 



which can be written in full as: 



2bw cos a 
—w + | ; v 



a + b 



cos a 



+ a/o 2 — sin 2 a) 



1 + a 



2 v a 2 — sin 2 a 2bw cos a 

> v. 

a a + b 



(A25) 

The factors (\/a 2 — sin 2 a) /a and (cos a + \/ d 2 — sin 2 a) / (1 + a) are less than or equal to 1 
so we can replace them with 1 and the inequality for the fifth condition remains true: 



2b ( 2bw cos a \ 2bw cos a 

■ v — w ] > v. 



(A26) 



1 + b \ a + b J a + b 

The following sequence of algebraic manipulations demonstrates the impossibility of the 
sequence (32313): 



2b ( 2bw cos a \ 2bw cos a 

— v — w ] > v 



1 + 6 V a + b 
2bw cos a (b — 1 



a+b V 6+ 1 



> 1 T+l ' " + 



a + b 
2bw 



1 + b 



2bw 

T+b 



(b — 1) cosa 
a + b 



> 1 FT! ' " > 



(b — 1) cos a > a + b 



(A27) 

(A28) 

(A29) 
(A30) 



which is a contradiction. 



Appendix B: Proof of Lemma 2 

Suppose particles 1 and 2 collide and that before the collision, the velocities of 2 and 3 
are p = p exp(i(j)) and q = q exp(iip). Then after the collision, their velocities, p' and q', 
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are given by: 

/ 




x \ I p \ I xp 

\(x-l)/(a + l) I J \q) \q+p(x-l)/(a + l) 



(Bl) 



After 1 and 2 collide, 1 and 3 must collide if we are considering a process of the form (32)™ 
or (2)(32) n . Now if the previous collision between 1 and 2 was preceded by one between 
1 and 3, 3 is moving radially away from 1 at the moment of the collision between 1 and 
2. So for another collision to occur between 1 and 3, the velocity of 3 after the collision, 
q' = q + p(x — l)/(a + 1), must be directed toward the origin and so be antiparallel to its 
velocity before the collision, q. But if q and q + p(x — l)/(a + 1) are antiparallel, the same 
is true of q and p(x — l)/(a + 1). The phase of a; is a and by a trigonometric identity, the 
phase of x — 1 is (a + 7r)/2. Therefore, for q and p{x — l)/(a + 1) to be antiparallel, we need: 

^ + vr = 0+(a + 7r)/2 (B2) 
a = 20 - 0) + 7T (B3) 

The above derivation was for the required scattering angle of a collision between 1 and 2 in 
order for a collision between 1 and 3 to follow. But since the final result is independent of 
any mass ratios, it also works for the case of a collision between 1 and 3. However for that 
case, by <fi we mean the angle of incidence of particle 3 and by ip we mean the initial angle 
of 2's velocity. 

Following the collision between 1 and 2 described above, the phase of 2's velocity is: 
■ip' = <p + a = 2ip — + 7r and the phase of 3's is = ip + 7r. Now if these are plugged into 
the formula for scattering angle, we obtain: 

a' = 2<y - 0') + 7T = 20 - 0) + 7T = a. (B4) 

This means that if four collisions occur in a row (either (2323) or (3232)) the second and 
third scattering angles are equal. The first and fourth collisions are required because the 
collisions discussed above were assumed to be preceded and followed by collisions. Processes 
which lead to capture leave 1 and 2 with the same velocity. Two particles with the same 
velocity may be said to collide after an infinite time so the process concludes with a virtual 
collision (3). Since this may be counted as the fourth collision, all collisions have the same 
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scattering angle with the exception of the first. 
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